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% Abstract 

The Higgs- like signal observed at the LHC could be due to several mass degenerate resonances. 
We show that the number of resonances is related to the rank of a "production and decay" matrix, 
Rif. Each entry in this matrix contains the observed rate in a particular production mode % and final 
state /. In the case of N non-interfering resonances, the rank of R is, at most, N. If interference 
plays a role, the maximum rank is N or in the CP limit N(N + l)/2. As an illustration we use 
the present experimental data to constrain the rank of the corresponding matrix. We estimate the 
^ LHC reach of probing two and three resonances under various speculations on future measurements 



and uncertainties. 
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I. INTRODUCTION 



At present there is clear evidence for a new boson at the Large Hadron Collider (LHC) 
with a mass of about 126 GeV [U |2]. In this paper we are interested in the possibility that 
the signal is due to two or more mass degenerate states. We are especially interested in the 
question whether one can tell experimentally if there are more resonances, and perhaps even 
constrain their number, without actually resolving the resonant peaks. That is, we are in- 
terested in the case that the resonances are mass degenerate to better than the experimental 
resolution. 

So far, the observed rates Rif of the 126 GeV resonance have been consistent with the 
Standard Model (SM) Higgs boson, although their central values do exhibit 0(1) devia- 
tions. Here i runs over the production modes: gluon fusion, vector boson fusion, associated 
production, while the final states / are 77, rr, etc. By itself, any deviation from the SM 
prediction in one element of Rif can be explained by some kind of a new physics scenario 
with one resonance. Some patterns of deviations in Rif, however, are bound to be incon- 
sistent with the hypothesis of a single resonance, even if one allows for non-SM couplings. 
This is easily understood using a counting argument. We denote the number of observed 
production (decay) modes by n, (n./). Assuming there is a single resonance, the unknowns 
are the rii cross sections {o"j} and rif branching ratios {Bf}. The number of observables is 
riirif, which are all the entires in Rif. Thus, the set of equations one has to solve is 

cnBf = Rif. (1) 

In this system of equations the number of unknowns, rii + rif, can be less than the number of 
observables, riirif. With enough measurements, this may lead to a system of equations which 
is overdetermined and potentially inconsistent [j] Such inconsistency cannot be removed by 
any type of new physics that affects only the production or decay of a single resonance. We 
would then be forced to abandon the single resonance hypothesis and conclude that there 
are more resonances. 

The idea of detecting two resonances using such inconsistency has been pointed out 
in pEJ [5] in the context of the NMSSM, and in [6] for two Higgs doublet models, where the 
authors discussed conditions that must hold in a single resonance scenario. In this paper we 
generalize these conditions to include any number of resonances. 

Assuming no interference between the resonances we find that rank(i?), the rank of the 
signal matrix Rif, constitutes a lower bound on the number of resonances required to explain 
the data. For interfering resonances we find, on the other hand, that this lower bound is 
relaxed to the integer that is greater or equal to y rank(i?) in the most general case or to 
(y/l + 8rank(i?) — l)/2 if CP is assumed. It follows that if rank(i?) > 2, there must be at 

1 For a single resonance, by taking the logarithm of each equation in ([I]) one arrives at a system of linear 
equations with logCi and logi?/ as the new variables. Therefore, discussing linear dependence, over- and 
undcr-dctcrmination, and (in)consistency is formally justified. 
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least two resonances, whereas if rank(i?) > 3, the resonances must interfere, or there must 
be at least three of them. 



II. GENERAL DISCUSSION 

We begin by recalling a few useful mathematical relations. Consider two sets of vectors 



in 



u r = (u[,...,u r N ), v r = (v[,...,v r N ), (2) 



where r = 1, . . . , n and we assume N ^> n. We define two N x N matrices whose elements 



are 

n 

|2 



A v = K v il ' B *i 



r=l 



E 

r=l 



r r 

u iVj 



(3) 



Then one can prove that generically, 

ranky4 = n, rank I? = n 2 . (4) 

When the phase between u\v r ^ and u\ is independent of i, j, the above results are modified 
such that 

n(n + 1) 

rankA = n, rank I? = . (5) 

Note that the condition leading to ^ is satisfied when the vectors are real. 

More generally, u and v can have different dimensions, and we may have cancellations if 
the elements are not generic. Thus a more general result is to treat the r.h.s. of the relations 
in Eqs. Q and ^ and as an upper bound on the corresponding ranks. Of course, the rank 
is also bounded by N. Equipped with these results we now return to physics. 

Suppose there exist n r > 1 resonances {</>i, • • • ,<fi nr } all with masses close enough so 
that the experimental resolution is not sufficient to resolve them. If the masses of different 
resonances are split by much more than their respective decay widths, interference effects 
are highly suppressed and can be neglected. The signal strengths Rif are then given by 

R lf = Y,°lB r f = *i-B f , (6) 

r=l 

where the production cross sections a\ = a^/af M and branching ratios B\ = Bl/Bf M for 
each resonance r are normalized to the SM values. The sum in Eq. ^ runs over the 
degenerate resonances r = 1, . . . , n r . Using Eq. Q we see that 

n r > rank(i?). (7) 

Namely, if one is able to measure rank(i?), it would imply a lower bound on the number 
of resonances. Recall that the rank of a matrix Rif is equal to the number of linearly 



3 



independent rows or columns, and may be obtained, e.g., using singular value decomposition 
(SVD), where the number of nonzero singular values equals rank(i2). 

If the decay widths are larger than the mass splittings between resonances, interference 
effects may be important. For instance, for perfectly overlapping resonances, 

2 

= \Ai-A s \ 2 , (8) 

where A\ and A r j are the properly normalized amplitudes for the production and decay of 
the resonance r, such that a\ = \A^\ 2 and Bj = \A T f\ 2 . In contrast to the non interfering 
case, the elements of Rif are no longer given by a sum of positive terms, and cancellations 
may occur. Thus, small entries in Rif are possible, even when a\ and B T f are large, providing 
more model building freedom. 

Moreover, interference affects the relation between n r and the rank of R. Indeed, applying 
Eq. Q to Eq. g results in 

rank(i?) < n 2 r n r > y / rank(i2). (9) 

In the case of CP conservation, the amplitudes for a given resonance have the same phase. 
This is true since the only significant CP even phase ("strong phase") is due to the decay 
width in the Breit-Wigner s-dependence of the decay amplitude and is universal. Non- 
universal effects such as rescattering phases are perturbative and can be neglected. Thus, 
in the CP conserving case we can apply Eq. (|5| to Eq. (18]) and obtain 

, / n r (n r + 1) . . 

rank(/r) < '-. (10) 

Solving for n r we obtain the lower bound 

-(y/l + 8rank(i?) - l). (11) 

It follows that any data which can be explained using n r non-interfering resonances may 
possibly also be explained with fewer interfering ones. For example, if the data satisfies 
rank(i?) = 3, there must be at least 2 interfering resonances, or 3 resonances without 
interference among them. In the case where the data suggest rank(i?) = 4, there must be 
at least 2 interfering resonances with CP violation, or 3 interfering resonances without CP 
violation or 4 non interfering ones. 

We conclude this section with a few remarks. 

1. First we comment on the transition between the interfering and non-interfering case. 
This transition must be smooth, yet our result above seems discontinuous. This may be 
understood as follows. When interference is taken into account, the extra non vanishing 
singular values of the matrix R are proportional to the size of the interference term. 



R 



if 



E 



a:a\ 
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prod.\decay 77 WW* tt ZZ* bb 

gg 1.6 ±0.35 PUCE] 0.8 ± 0.3 [13 EH] 1.2 ± 0.8 £E21 E3] 1.0 ± 0.3 PS] 

VBF 2.1 ±0.9 [mE] -0.2 ± 0.6 [13 EH] 0.3 ± 0.7 [m H3] 

FiJ 1.9 ±2.6 pi] -0.3 ± 2.1 jH] 1.0 ± 1.8 [S] - 0.8 ± 0.6 [H] 

ftfc - - < 3.8 [IS] 

TABLE I. The experimental values of signal strengths, Rif, for several production modes (listed 
in the first column) and different final states (listed in the first row) obtained by naively averaging 
ATLAS and CMS results with statistical and systematic errors added in quadrature and neglecting 
correlations. 



When the resonances are far apart, their interference is negligible and proportional 
to their overlap, and the singular values approach zero. For example, consider two 
resonances with equal width T and mass difference Am. In this case, we expect two 
singular values to be similar, while the third one is suppressed by the overlap. In the 
Breit-Wigner approximation, it is given by (2r/Am) 2 , although realistic profiles may 
deviate from the Breit-Wigner curve at their tails. 

2. In the limit of conserved CP the interference can still increase the rank of R, however, 
the resulting matrix is not generic. For example, when all the amplitudes are real, the 
matrix = \J Rij has the same rank as that of Rij in the non interfering case. 

3. In all of the above, we have ignored effects of interference between the signal and the 
SM background. Such effects were studied, for example, in the case of gg — > H — >■ 
77 [3 IS]- While this effect can be experimentally taken into account in the SM, it 
cannot be done model independently. Yet, taking the SM as a benchmark, we may 
estimate the effect to be small enough, and ignore it. It may become an issue in the 
future. 



III. NUMERICAL EXAMPLES 



We discuss next a few concrete numerical examples, starting with the information that can 
be extracted already with present experimental precision and then contemplate projections 
to the future. The current status of the signal strength matrix Rif is given in Table |TJ where 
we have averaged the experimental results from ATLAS and CMS. 

In general, the signal strengths Rif are measured with finite precision 8 Rif. We simplify 
the discussion by neglecting correlations between different channels, cross contaminations 
and other systematic errors. For instance, depending on the cuts, some of the gg signal 
can appear in the vector boson fusion (VBF) production channel, as is the case in h — > 77. 
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These effects were partially taken into account in Table [IJ e.g. for the 77 final state. We also 
assume that the acceptances in the NP models are similar to those of a SM Higgs boson. 
This may be simply realized by assuming that different resonances are described by the 
same low energy effective Lagrangian as the SM Higgs J2HH23]- While these approximations 
suffice at present precision, we emphasize that these effects can be important and there is 
no reason for them not to be included in a full-fledged analysis. 

To constrain the number of resonances, we are interested in the statistical properties of 
the hypothesis that rank(i?) > r. Note that one can only put a lower bound on the rank 
since including experimental errors generally increase it. Finding the significances for each 
possible r may be done in a number of ways. Below we discuss two of these, but there 
are many other possibilities. The best approach then depends on the particular situation. 
Depending on the exact situation, one way may be better than another. 



A. Using Singular Value Decomposition 

Obtaining the rank of R may be done most generally by using singular value decompo- 
sition (SVD), 

R = USV T } (12) 

where U and V are orthogonal matrices with dimensions n,i and rif respectively, and S 
is a rectangular (n» x nf) diagonal matrix. The diagonal elements of S are the singular 
values, Aj = Su (no summation), and may be obtained by diagonalizing RR T (using U) or 
R T R (using V). By definition Aj are non-negative, which can always be accomplished by 
appropriate redefinition of U and V. It is convenient to define S, U and V, such that the 
singular values are ordered with Ai the largest. The rank of an exactly known matrix is 
given by the number of nonzero singular values N(Xi ^ 0). Since in our case the matrix R 
is only known to a certain precision, these singular values have distributions induced by the 
uncertainties of the matrix elements. To obtain a statistical test of rank(i?) we follow the 
procedure of [24]. We split U, S and V into the block matrices 





$= 1 ... 1 . v= \ Vr " : 1 . us) 



' rxr 

V 



where U is an (n» — r) x [m — r) square matrix, S an (n» — r) x (nf—r) rectangular diagonal 
matrix, and V and (nf — r) x (nf — r) square matrix. We then construct 

A r = USV T } (14) 

which is an (n, — r) x (nf — r) matrix that is zero if the singular values are zero. (Note that 
we differ in the normalization of A r from [21] •) We can then test whether singular values 
are consistent with zero simply by asking whether A r is consistent with zero. We neglect 
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A r mat. el. val. sing. val. 



FIG. 1. Left: The probability distributions for the (Ai)n (blue) and (Ai)i2 (red) entries with 
present data and for 10 x more data assuming present central values and the dominance of statistical 
errors. The Rif matrix has entries from i = {gg,V15F} and / = {77, rr, WW} channels. Right: 
The probability distributions for the two singular values of the 2 x 3 Rif matrix. 



correlations and assume that the errors for the matrix elements of A r are Gaussian, which 
we denote by cr[(A r )jj]. We then define 

which approaches the x 2 distribution for (rii — r) x {rif —r) degrees of freedom in the large 
statistics limit. While we do not propagate the correlations, one could easily modify the 
above definition of x 2 (A r ) to take into account the correlations. On average we expect 
X 2 (A r ) = (rii — r) x (rif — r). Much higher values would correspond to A r significantly 
nonzero, with CL obtained from x 2 distribution, signaling that the rank of the matrix R is 
at least r + 1. The advantage of working with A r instead of with singular values is that the 
entries of A r are not bound to be non-negative (cf. Fig. [T]) so that it is easier to devise such 
a simple x 2 t es t- 

Next, we apply the above method to the data. The largest matrix Rif that we can 
form using present measurements is 3 x 3, with i = {gg,VBF,VH} and / = {77, rr, WW} 
collected in Table HJ Since the last row is least well measured we choose to use the 2x3 
rectangular submatrix i = {gg,VBF} and / = {77, rr, PFPF} for the numerical examples 
(we will return to the 3x3 case below). The most obvious test is for rank(i?) > 1, choosing 
r = in (13). This means that A = R. Neglecting correlations we have x 2 (Ao) = 36.1 for 
6 d.o.f., which corresponds to a probability 1 — CL = 3 ■ 1CT 6 that this is just a statistical 
fluctuation. The signal significance is 4.7a and we confirm the hypothesis of at least one 
resonance. 

Testing whether rank(i?) = 2 we take r = 1, so that A x is a 1 x 2 matrix. The probability 
distributions for the two entries in Ai are shown in Fig. [I] (left) and were obtained from simple 
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Monte Carlo with 10 6 iterations assuming Gaussian errors for the Rif matrix elements. This 
gives x 2 (Ai) — 0.8 which is not statistically significant for 2 d.o.f. as it corresponds to 
0.4cx signal significance. That is, with the current data we cannot test for more than one 
resonance. With ten times more data, however, assuming the same central values and scaling 
of errors as 1/y/L, one would obtain a signal for rank(i?) = 2 with 3.9a significance. While 
this is just a very crude estimate, we conclude that with roughly an order of magnitude 
more data we will be able to probe for a rank two matrix. 



B. Using The Determinant 



While the method described above is very general, in some cases we can devise simpler 
tests. When R is a square block matrix (rij = nf — n) one can check if det R is significantly 
nonzero, which would imply that the rank of R is n. For 2x2 blocks (i.e., rii — nf — 
2), the determinant already contains the maximal information available on the number of 
resonances, since det R = implies that one resonance suffices, while a det R 7^ implies 
n r > 2. For larger square blocks, however, the determinant can only distinguish the n r > n 
case against n r < n. We will now give examples for both 2x2 and 3x3 blocks. While the 
former probes two resonances, the latter allows to probe three resonances and interference. 

First, consider the two 2x2 sub-matrices R = ({gg, VBF}, {77, rr}) and R B = 
({gg> VBF}, {77, WW}). The SM is excluded if any R±j 7^ 1. The single resonance hy- 
pothesis is excluded if det R A 7^ or det R B 7^ 0. When formulated in term of the ratios 

RvBF,WW Rgg,WW -RvBF,rr ^gg,TT M c n 

n = -5 , r 2 = , r 3 = , r 4 = , (16) 

-'•'VBF, 77 99,77 -**'VBF,77 99,77 

these inequalities become r\ 7^ r 2 or r 3 7^ r^. Ratios such as above have been considered 
in [1H6], and have the virtue that certain systematic errors cancel [9]. As an example, 
consider a situation in which all the rates remain at their current central values with much 
reduced errors, except -Rvbf ,wwi which we take to be small instead of the current negative 
value: 

"^99,77 = 1-6) -RvBF,77 = 2.1, Rgg,WW = 0.8, RvBF ,WW = 0.1. (17) 

Interpreting such hypothetical results as coming from a single on-shell particle is inconsistent, 
since the resulting signal matrix has nonzero determinant, detR 7^ 0. In terms of ratios, 
this fact is evident in the inequality of the ratios 

R ™> ww = 0.048, %^ = 0.5, (18) 

-'•'VBF ,77 99 ,77 

while for a single resonance the two ratios would be identical and equal to Brwwl -Br 77 . If 
this inequality becomes certain in the future, it will therefore point to rank(i?) = 2, implying 
the existence of at least two mass degenerate resonances, 4>i an d <t>2- Such hypothetical result 
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0.0 0.5 1.0 1.5 2.0 -0.5 0.0 0.5 1.0 1.5 2.0 

fl—R VBF, Ww/RvBF,yy >~3 =RvBF, tt/R VBF.yy 



FIG. 2. Contours of constant integrated luminosity at LHC 14 TeV required to exclude at 3<r 
the SM (the single resonance hypothesis) if the measurement is within blue (green) regions in the 
(ri,r2) plane (left) and in the {r^^r^) plane (right). We have assumed statistical scaling of errors, 
while present la measurements are denoted with black error bars. The 14 TeV run projections 
were obtained by rescaling the current 7 + 8 TeV errors with the increase in the higgs production 
cross section [25] to reflect the gain in statistics. 



could be explained, for example if 0i does not couple to W bosons while 2 does not couple to 
photons. Note that since we have assumed low R\bf,wwi if the resonances do not interfere, 
some of the parameters need to be suppressed, for example, the "Higgs" coupling to WW. 
This is because in that case the Rif values are sums of positive numbers. In contrast, no 
such suppression is needed for interfering resonances where the suppression can come from 
negative interference. This provides more model building freedom. 

A more realistic discussion includes experimental errors, where one tests for a statistically 
significant deviation from det R = 0. In Fig.[2j we show contours of the integrated luminosity 
required to exclude at 3a the SM (blue) or the single resonance hypothesis (green), in the 
(rq, r 2 ) and {r^^r^) planes. In the extrapolation we assume statistical scaling of errors, that 
is oc 1/ y/L, and central values taken from Table [ij We find that if the current central values 
remain where they are now, then O(300fb _1 ) at LHC 14TeV would suffice to exclude at 3<r 
the single resonance hypothesis. 

Finally, we discuss the possibility of testing for rank 3. Note that this is not only a 
generalization from 2 resonances to 3, it is also the lowest possible rank which can probe 
interference. As an example we test for rank 3 using det R ^ test for n = 3. This distin- 
guishes between rank(i?) = 3 and rank(i?) = 1,2 cases. The rank(i?) = 1 and rank(i?) = 2 
cases may be resolved using SVD, as discussed in Eqs. ( 12 )-( 15 ). For inputs we use the 
({gg, VBF, VH}, {77, WW, tt}) block in Table |l| 
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First, consider an idealized future scenario, where all the central values are as in Tabled], 
except the negative ones, which we take to be +0.1, for illustration purposes. The relevant 
matrix has rank 3, and can be decomposed as usual using three non-overlapping resonances. 
The same matrix can also result from 2 resonances, provided they interfere (cf. Sec. [n]). For 
example, the two resonances may have the following amplitudes, 

Agg = (1, 0.458 e°- 5i ) , A VB f = (0.673 e" 1 ^, 1.088 e 0Al ) , A VH = (l.ll e 05i , 1.277 e 2i ) , 
A 77 = (l,1.95e L2i ) , A ww * = (0.9 e L9 \ 0.69 e" 314 ) , A Tf = (l.l e 0Ai , 0.17 e L6i ) . (19) 

Note that this example employs CP violation. As we mentioned before, a CP conserving 
example does not always exist. However, many times a CP conserving model can produce a 
matrix that is in the vicinity of the one that is needed. Thus, it is not clear how much data 
one will need in order to establish CP violation. 

We are now ready to discuss the more realistic case with finite statistics. We show the 
projected sensitivity in Fig. ^/e/i), assuming statistical scaling of errors, and assuming the 
cross sections scale with center of mass energy as they do in the SM. In this case C(300fb _1 ) 
at LHC 14 TeV are needed for 3a evidence of two resonance case using the rank(i?) = 2 test 
described above. We estimate that C(2000fb" 1 ) at LHC 14 TeV are needed for 3a evidence 
of three overlapping resonances using the test that deti? ^ 0. It is important to note 
that all these estimates crucially depend on the input parameters and 0(1) change in the 
central value of just one parameter can lead to a factor of few smaller required luminosities. 
Therefore, in the right panel we show a plot where the central values are allowed to vary. 
Contours of constant luminosity are shown where the significance for discovering 3 non- 
interfering or 2 interfering resonances (det R ^ 0) reaches 3a. The error has been propagated 
linearly by expanding the determinant. While this is clearly not adequate for large relative 
errors, it provides a rough estimate for the significance. 



IV. SUMMARY AND OUTLOOK 

We have shown that one can test for the presence of multiple resonances in the recently 
observed Higgs signal, even if the multiple resonances are degenerate in mass with respect to 
the experimental resolution. The test requires measuring algebraic properties of the event 
rate matrix, R, in particular its rank. 

For purposes of the current Higgs signal at the LHC, we find that if the current central 
values remain then O(100fb _1 ) - O(300fb _1 ) of integrated luminosity at the 14 TeV LHC 
may be sufficient to detect evidence for 2 resonances with 3a significance. Detecting three 
resonances (or interference between two) would require more data, C(2000fb~ 1 ). Of course, 
the situation can dramatically change based on how the central values move. 

While in this work we pointed out that the rank of the R is an important observable, we 
did not carry the full statistical treatment of the available data. This task must be performed 
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FIG. 3. Left: The projected significance for a discovery of 3 (2) overlapping resonances using 3x3 
R matrix (see text) using the det R / (rank(i?) > 2) tests are shown with solid and dashed blue 
(red) lines for LHC at 8TeV and LHC at 14 TeV, respectively, as a function of integrated luminosity. 
Right: Integrated luminosity at 14 TeV LHC required to probe 3 non-interfering or 2 interfering 
resonances with 3<r significance as a function of the ratios R g9: ww* I 'Rgg,rf and Rvh,ww* /Rvh,tt> 
keeping R gg> ww* an d Rvh,tt fixed to their central values in Table 11} 



by the experimental collaborations, taking into account correlations and systematics, and 
correctly accounting for the propagation of error. 

The principles we have discussed in this work are universal and may prove suitable not 
just for Higgs physics at the LHC. For example, new resonance peaks may be found in the 
following years, as motivated by many extensions of the SM, and investigating their signal 
matrices might be the only way to probe degeneracy, especially in a hadron collider where 
the mass resolution cannot be very good. 
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